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Using Monte Carlo simulations, we have studied aging phenomena in three-dimensional 
Gaussian Ising spin-glass model focusing on quasi-equilibrium behavior of the spin auto- 
correlation functions. Weak violation of the time translational invariance in the quasi- 
equilibrium regime is analyzed in terms of effective stiffness for droplet excitations in the 
presence of domain walls. The simulated results in not only isothermal but also T-shift 
aging processes exhibit the expected scaling behavior with respect to the characteristic 
length scales associated with droplet excitations and domain walls in spite of the fact 
that the growth law for these length scales still shows a pre-asymptotic behavior com- 
pared with the asymptotic form proposed by the droplet theory. Implications of our 
simulational results are also discussed in relation to experimental observations. 
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1. Introduction 

In recent years aging dynamics in spin glasses has been extensively studied. BH a 
basic experimental protocol is isothermal aging process which is relaxation after the 
spin-glass system is quenched from temperature above the spin-glass temperature T c 
to temperature T below T c . One of the most interesting results in the simulational 
studies on 3D short-ranged Ising spin-glass model is that the coherence length Rr{t) 
of the correlation function of real two replicas at time t after the quench grows in a 



power-law as,H'En3 (but see alsoQ) 

R T {t)~L {t/T ) l ' z , (1.1) 

where Lq is a certain characteristic length and exponent 1/z depends an T. The 
latter was obtained as l/z(T) ~ 0.17T for T < 0.7 in our previous paperQ which we 
refer to I hereafter. It is tempting to regard this Rr(t) as a characteristic length 
of the ordering process, i.e., a mean distance of domain walls separating different 
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pure states, which are parallel or anti-parallel to a ground state. In this respect, it 
is only the droplet theorytm at the moment, that provides us some concrete scaling 
arguments based on such characteristic length scales as Rrit). 

According to the droplet theory by Fisher and Huse (FH),Q isothermal aging 
in spin glasses is associated with coarsening of domain walls, which is driven by 
successive nucleation of thermally activated droplets. Up to waiting time t w after 
the quench, domains of averaged size -Rt(^w) have grown up. Within each domain 
small droplets of size L (<C -Rt(^w)) are thermally fluctuating as they are almost in 
equilibrium: a typical value of their excitation gap F^ p scales as 

Ff p ~ T(L/L Q ) e , (1.2) 

where T is the stiffness constant, and that of free-energy barrier B^ p scales as 

~ A(L/L f, (1.3) 

where A is a characteristic free-energy scale. The above two exponents satisfy 
9 < (d—l)/2 and 9 < ip < d — 1 with d being dimension of the system. 

However, there are some droplets which touch with the domain wall so that their 



excitation gap is reduced from eq.(1.2). This effect was studied by FH by making 
use of some geometrical and probabilistic arguments. Averaged over all droplets 
including those within the domain, FH has derived a typical value of excitation gaps 
of droplets with size L in the presence of the domain walls with mean separation R 
as 

F typ R = T ^ [L/R]{L/Lo) e : (L4) 

with the effective stiffness constant given by 

T cS [L/R] = T^l-c v ^ d \ (1.5) 

where c v is a numerical constant. 

By making use of the above argument on the effective stiffness constant, the 
spin auto-correlation function 



C T (T;t w ) = Ci,T(T;t v ) with Ci,r(r;i w ) - (Si(r + t w )S i (t w )) T , (1.6) 

is evaluated as the following. Some further details of its derivation are described in 
our separated papenlJ which we refer to II hereafter. A droplet of size L enclosing 
spin Si contributes to C^rCntw) as 

r fw ^ ~ / 1 for TL ^ > r ' n 7^ 

l ' T[ ' wj " 1 <Si> 2 T ~ 1 - 4exp(-F L , fl ( i )/T) for T L (i) < r, ^ 

where Tl(i) = roexp(_Bi(i)/T) is the relaxation time of the droplet whose free- 
energy barrier is Bl(i). The averages (denoted by the over-line in eq.( |l.6| )) over 
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site i and different realizations of interactions (samples) are all taken into account 
by the probability distribution of free-energy gaps F^^i) which is assumed to scale 
as. 



Pl.r(F) 



1 



F 



(1.8) 



L,R \ ± L,R/ 

where a scaling function p(x) satisfies p(0) > 0. Then, gathering the multiplicative 
contributions of droplets with various sizes, L — 2 ,1 Lq; n — 0, 1, 2, enclosing spin 
Si, we obtain 



l-C T (r;t w ) 



Lt(t) dL p(0)T 



1 



L (L/L )° r° s {L/R T (t w )} 



1 - C cq , T (r) + 



Lt(t) 

T(L T (T)/L ) e \R T (t v ) 



Cl/ 5(0)T 



d-e 



where c\ is a numerical constant and 

Tp{0) [ Lt[t) dL I 



1 - C cq , T (r) 



Here Rr(t) is given by 



Tm 

L (L/L ) 9 6T 



R T (t) ~ L [(r/A)ln(i/ro)] 1/V; . 



L 



L t {t) 



(1.9) 
,(1.10) 

(1.11) 
(1.12) 



and Lt(t) by the same expression with t replaced by r in accordance with eq.( |1.3| ). 
The latter is the characteristic length of droplets which can fluctuate within time 



scale t at temperature T. The above scaling expression eq.(l.K) is appropriate in 
the regime Lt(t) <C i?r(^w) which defines the quasi-equilibrium regime of present 



interest. The third term in eq.(l.lO) represents the leading correction to equilibrium 
behaviour which weakly violates the time translational invariance in this regime. 



The simulational result eq.(l.l) differs from eq.(1.12). A possible interpretation 
of this discrepancy which we will pursue in the present work is as follow. The 
difference in the growth laws of Rrit) and Lt(t) is simply due to two different 
time windows which one observes the aging processes; one is an asymptotic regime 
for which the droplet theory is constructed, and the other a pre-asymptotic regime 
that the simulations can observe. However, we consider that the scaling expres- 



sions written in terms of Rr(t) and Lt(t), such as eq.( l.lC ), are common to both 
regimes. A similar interpretation has been already proposed for some aspects of 
aging dynamics.00 In our paper II we have shown by detailed analysis on the spin 
auto-correlation function that the droplet picture in the sense described above in 
fact holds in the quasi-equilibrium of the isothermal aging in the 3D short-ranged 
Ising spin-glass model. 

The main purpose of the present work is to demonstrate that the same pic- 
ture consistently describes behavior of the aging dynamics in the quasi-equilibrium 
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regime of the so-called T-shift process. i'El The latter is another typical experimental 
protocol of aging: first to let the system age isothermally at T\ , to change temper- 
ature to T2 at waiting time t w i , to let the system age isothermally at T2 in a period 
t W 2 at which a small probing field is added, and then to measure response of the 
system at time r afterwards. We may write the corresponding correlation function 
as CT 2 ,Ti(i~;t W 2;i w i)- 

One of most interesting problems in the T-shift process is whether or how we 
can observe chaotic nature of the spin-glass phase which is derived from the droplet 
theory;ElEj namely, there exists such a characteristic length scale Li^-t^, called the 
overlap length, that the equilibrium spin-glass orders at T\ and T 2 differ entirely 
from each other in the length-scale larger than L\t 1 -t 2 \- I n m e present work, as for 
the first step of our study on this problem, we have analyzed Ct 2 ,Ti (t; i W 2 ; twi ) with 
tw2 = in the quasi- equilibrium regime which is specified by r <C t w \. Up to now 
even a precursor of the chaotic nature has been ascertained within the time scale of 
our simulation. Instead, as mentioned above, the simulated results are consistently 
interpreted by the droplet picture without introduction of L\t 1 -To\- I n other words, 
L|T 1 _y 2 |, if it exists, is much larger than the length scale of our simulation. 

After describing our model and numerical method in the next section, we briefly 
review the results on behaviors of the correlation function in the quasi-equilibrium 
regime of the isothermal aging in §3. We then present and discuss the corresponding 
results in the T-shift aging process in §4. In the last section we argue implications of 
our simulational results on the experimental observations based on the fluctuation- 
dissipation theorem. 



2. Model and Method 

We have carried out standard heat-bath Monte Carlo simulation on aging phe- 
nomena in the 3D fsing spin-glass model with Gaussian nearest-neighbor interac- 
tions with zero mean and variance J = 1. The spin-glass transition temperature 
is numerically determined most recently as T c = 0.95 ± 0.04.E-3 The data we will 
discuss below are obtained at T = 0.5 ~ 0.8 in L s — 24 systernaveraged over 160 
samples with one MC run for each sample. In our previous work,til hereafter referred 
to I, it was confirmed that finite-size effects do not appear within our time window 
( < 2 x 10 5 MCS) for these parameters. Also in I we showed from the fmite-size- 
scaling analysis of relaxation of the energy per spin in isothermal aging that the 



exponent 9 in eq.(1.2) at T = 0.7, 0.8 is estimated as 9 = 0.20±0.03 which coincides 
with the result of the defect energy analysis at T = 0.E3 



3. Results on Isothermal Aging 

It is now well established that the correlation function Ct(t;£ w ) in isothermal 
aging consists of two characteristic time regimes, one is the range t «i w which 
is the quasi-equilibrium regime mentioned in §1, and the other r 3> t w called as 
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the out-of-equilibrium (or aging) regime. Our interest here is in its behaviors in 
the former regime. As shown in our paper II, characteristic features of Ct(t;<: w ) 
in the quasi-equilibrium regime are more easily understood when 1 — Ct(t;£ w ) 
are analyzed as a function of t w /r with r being fixed. It turns out that such data 
points of different r lie top on each others when they are shifted vertically by proper 
amounts qt(t). The results of such a one-parameter scaling are demonstrated in 
Fig. 1 where we plot 

AC T {T;t w ) = l-C T (T;t w )-a T (T) (3.1) 

against t w /r. The {ay(r}} thus determined are presented in the inset of Fig. 1 
where we also show the results of the same analysis at different T's. They appear 
as almost linear functions of lnr. 



In order to compare the above results with the scaling expressions eqs.( l.lC ) and 



( |f.ll|) , we only need to notice the fact that factor (L /Lt(t)) can be practically 



expanded as 



1 ' M 1 - -^ln - , (3.2) 



L t (t) ) \T ) z{T) \t 



since 9/z(T) is small as compared with l/ln(r max /ro), where r max (— 512) is the 
maximum t in the present observation where we put tq = 1. Indeed, 9/ z(T) ~ 0.02 
with 9 = 0.20 ± 0.03 and l/z(T = 0.6) ~ 0.102, while l/ln(r max /T ) ~ 0.160. This 
factor in cq.( l.lOj ) can be approximated by its leading term (= 1), so that the last 



term in eq.(LK) becomes a function of only Lt(t) / 'it/r(£w) , which is given by 

AC T (r;t w ) oc {L t {t) / R T (t w )) K , (3.3) 
with k = d — 9. This is consistent with the fact that the one-parameter scaling 



shown in Fig. 1 does work well. The leading term of eq.(l.ll), on the other hand, 



comes out from the second term in eq.( p.2| ), and we obtain 

a T (r) = 1 - C oq , T (r) ~ ^M_L_ i n ( L ) . , :U , 



T z{T) \t q 

This is quite consistent with nearly linear increase of «t(t) shown in the inset of 
Fig. 1. Furthermore we obtain daT{T) / '<91nr oc T 2A from the data at four tem- 



peratures shown in the figure. This is compatible with eq.(3.4) which predicts the 
corresponding slope proportional to T 2 (note that l/z(T) oc T). The extra In- 
dependence, which yields the power T 2A larger than T 2 , may be attributed to that 
of the factor p(0)/T. 
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Fig. 1. The function ACt(t;( w ) as a function of t w /r obtained in isothermal aging at T = 0.6. 
In the inset we show {a(r)} thus extracted at T = 0.5, 0.6, 0.7 and 0.8 from bottom to top. Four 
symbols at T = 0.6 represents four sets of {a(r)} which yield significantly different values of k in 
eq.@. 



As for the determination of value k in eq.(3.3), however, a more careful treat- 
ment is required since our simulational time window is not enough wide to extract 
precisely the limiting behavior of ACt(t;£ w ) at Lt{t)/ '-Rr(iw) — > 0. By the de- 
tailed analysis described in II, the values of k acceptable within accuracy of the 
present simulated data arc estimated as n ~ 2.3 ~ 3.1. The result is compatible 
with the expected value d — 9 ~ 2.8, though the data of ACr(T;t w ) alone cannot 
specify the value of 9. Let us remark also that this ambiguity in n little affects the 
values of {o:t(t)}, in particular, their slopes with respect to lnr. This is demon- 
strated in the inset of Fig. 1 by drawing four sets of {oct(t)} at T = 0.6 which yield 
K ~ 2.2,2.6,3.0, and 3.3, respectively. These circumstances are similar at other 
temperatures we have examined. Therefore we show the representative results of 
{ar(r)} in the inset of Fig. 1 which give rise to n ~ 2.8 ~ 3.0. 

So far we have demonstrated that our simulated data of Ct(t;^ w ) are consis- 



tently interpreted by the scaling ansatz eq.(1.10), when the growth law of Rr{t 



eq.(1.12) in the droplet theory is replaced by eq.(l.l). We regard that the former 
describes the growth law in the asymptotic regime closer to equilibrium, while the 
latter does in the pre-asymptotic regime, and that the common scaling forms such 
as eq.(l.lC) hold in both regimes. 



4. Results on T-Shift Process 

Among the T-shift aging processes described in §1, here we restrict ourselves to 
analysis on Ct 2 , T t (t; t W 2', twi) with t W 2 = 0, which we abbreviate as Ct^,Ti{t; i w i). 
At time domains with mean separation i?Ti(^wi) have grown up. It has been 
confirniedta that, by the T-shift to T2 at t w i, the correlation length of the replica- 
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overlap function does not decrease. Rather it continues to increase though relatively 
gradually (rapidly) when T\ > (<) T2, and finally it merges to the isothermal curve 
Rt 2 (t) at T2 at around r = tf^ which we define by the condition 

,z(T 2 )/*(T0 



Rr 2 (t 



efi n 

wl / 



at,, 



(4.1) 



with a being a numerical constant rather close to unity. This means that i?Ti(^wi) 
with t w i we have examined does not exceed the overlap length Li Tl _ T i mentioned 
in §1. That Lyr^-T^U m case it exists, is much larger than the length scale of 
the presently available simulation agrees with the previous workB In the following, 
therefore, we analyze our data discarding the chaotic effect. 

Then from the observed results mentioned above, we can simply think of a 
scenario that, in the time range r <C after the T-shift, there exist domains 
nearly in equilibrium at T\ and with mean separation of Rt ± (t w i ) , and droplets 
up to mean size of Lt 2 (t) are fluctuating by thermal noises of temperature T%. In 
this situation with i?Ti(^wi) 3> Lt 2 {t), which we call the quasi- equilibrium regime 
in the present T-shift aging process, -Rri(iwi) plays a role of i?r(iw) in equations 
which correspond to eqs.(l.lO), (1.1C) and (3.3). More explicitly, we expect that, if 
1 — Ct 2 ,Ti (t; t w i) for various r are plotted against t^/r, they look quite similar to 
the 1 — Ct(t;£ w ) versus t w /r. Furthermore, as is the case for the latter shown in 
Fig. 1,1 — Ct 2 .Ti (t; t w i) are expected to lie on a universal curve by vertical shift, 
0^2,21(7"), of the data for each r. Indeed this is the case as shown in Fig. 2 where 
we plot 

AC T2 , Tl (r;twi) = [1- CT 2 ,TAr;t wl )} - a T2 ^) (4-2) 

against t^/r. 




10000 



Fig. 2. The plot ACt 2 ,t x (t; t„i) versus t^/r for the T-shift process with Ti = 0.8 and T 2 = 0.6. 
The functions ft;; (t) thus extracted are presented in the inset: symbols with the dotted line 
represent ftT 2 =0.6,T 1 =o.8( r ) an d those with the broken line ftT 2 =0.8,T 1 =o.6( T )i while only the 
symbols do the isothermal arCO with T = 0.6, 0.7 and 0.8 from bottom to top. 
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The vertical shifts aT 2 .T 1 {'r) in eq.(4.2) depend on both T\ and To ex 



plain this let us go back to an equation for the local correlation function, which 



corresponds to eq.(1.7) for the isothermal process. In the present T-shift process 
deviation of the function from unity occurs also when a droplet enclosing spin Si 
has relaxation time tl(i) which is smaller than r, but the value it relaxes to is now 
given by 

QiTx.iifctwi) -< Si > Tl < St >t 2 - l-2[exp(-F L . fl (z)/T 1 ) + exp(-F L . fl (z)/T 2 )], 

(4.3) 

for t > tl (i) , since Si here is regarded in equilibrium of temperatures T\ and T 2 at 



times r = and r > Tz,(«), respectively. Therefore we obtain, instead of eq.(3.4), 



This expression, combined with eq.( fc.4j ), consistently explains the relative orders 
of slopes da(r) I dhiT (= Aa) of various a(r) shown in the inset of Fig. 2; namely, 
Aa .s > Aao.8,0.6 > Aa .7 > Aa .6,o.8 > Aa . 6 - 




0.001 1 ■ ' ' ■ ■ ■ <- 

Lrl^/RrlKM-fyRAK) 



Fig. 3. The double logarithmic plot of AC't 2 ,t 1 (t; iwl) versus Lt 2 { t )/Rti (twl) of the T-shift 
process for different r's are plotted by the different symbols (with Ti = 0.8 and T2 = 0.6). For 
comparison ACV(t; i w ) versus Ly(r)/i?j'(t w ) of the isothermal aging for different r's are plotted 
by the different curves (at T = 0.6). The slope of the dotted line is 3.0. 

Finally we show in Fig. 3 the double logarithmic plot of ACt 2 ,Ti ( T i *wi) versus 
(r/t^) 1 ' 2 ^ 2 ' which is proportional to Lt^^/R^ (i w i)- For comparison, we also 
plot ACt(t;£ w ) versus (r /t^ 1 ^^ (oc Lt{t)/ Rrit^)) of the isothermal aging. 
Similarly to the latter case, there remains a relatively large ambiguity in the value 
of the slopes, i.e., the exponent which corresponds to k in eq.( |3.3| ). But the result 
is also compatible with k = d — 8 within accuracy of the present simulations. 
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5. Discussions 

Our simulational results on the correlation function C(t,t') (— Ct(t;^ w ) with 
t = t + t w and t' = t w ) so far described are closely related to the response function 
G(t,t') via the fluctuation-dissipation theorem (FDT) 

GM^^Ml, (5.1, 

and so the magnetic susceptibilities frequently measured in experiments. That the 
above FDT holds at least approximately in the quasi-equilibrium regime even with 
weak violation of the time-translational invariance (TTI) has been carefully checked 
experimentajjy.0 It has been also demonstrated in the simulational work by Franz 
and Ricger,E£l though the authors did not mention it explicitly: their data which 
exhibit the FDT violate the TTI, i.e., depend also on t w . 



From eq.(5T) the following relation between the out-of-phase component of the 



ac susceptibility Xj-(w;t w ) and Ct(t w ;£ w ) in the isothermal aging is derived, 

(5.2) 



7T 3 



where r w = 2tt/lo. Its derivation and implication on the experimental results, such 
as the wi w -scaling of x"t(w; t w ), are described in II. 



Here let us emphasize that from eq.(5.1) we obtain also the relation between 
C(r;t w ) and the zero-field cooled susceptibility X ZFC ( T ;*w) as 

x ZFC ( T ;t w )= +W dt / G(r + t w ,t / ) = ^[l-C(T;t w )}. (5.3) 

Thus 1 — Oy(T;i w ) in the isothermal aging, which we have analyzed in §3, is di- 
rectly related to Xt ( T !^w) m the quasi-equilibrium regime. The similar relation 
is expected to hold in the quasi-equilibrium regime of the T-shift process which 
we have discussed in §4, i.e., when the response is examined by adding a probing 
field from t W 2 = 0. Although ~x^ c Jj]t w ) has been measured frequently since the 
pioneering work by Lundgren et f?|p and Xt^Ti ( t ! ^w) m the T-shift process with 
tw2 — has been also examinedJ13 most of the measurements have been focused 
on its crossover behaviour between the quasi-equilibrium and aging regimes. We 
believe that the detailed analyses on these x ZFC ( r ; *w) are of quite importance since 
they will present us evidences for (or against) the droplet picture we have discussed 
in the present work. 

To conclude, we have simulated aging dynamics in the 3D Ising spin-glass model, 
and analyzed behaviors of the spin auto-correlation function in the quasi-equilibrium 
regime of the T-shift aging processes. The simulated results, which are considered 
to represent aging properties in the pre-asymptotic regime, do not exhibit even a 
precursor of the chaotic nature of the spin-glass phase. Instead, they are satisfac- 
torily explained by the scaling argument based on the characteristic length scales 
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-Rri(twi) and Lt 2 (t), i.e., that of domains grown up till t w \ at T\ and that of 
droplet excitations in time scales of r at T2- However their growth law itself differs 
from the prediction by the droplet theory which is constructed for the asymptotic 
regime. We have conjectured that the experimental observation on X ZFC (t;£ w ) in 
the quasi-equilibrium regime is of quite interest to check this droplet picture. 
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